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I . Introduction 


Knowledge  of  the  blast  field  produced  by  Army  guns  is  becoming 
increasingly  important  as  attempts  are  made  to  increase  range  and  muz- 
zle velocity.  While  these  attributes  are  important  to  the  lethality 
of  the  weapon  system,  the  accompanying  side  effect  — stronger  muzzle 
blast  field  — can  prove  detrimental  to  the  launch  platform  as  well  as 
the  firing  personnel. 

Heretofore  information  on  the  blast  field  strength  was  obtained 
from  either  actual  firings  (expensive) , scaling  parameters  (always 
suspect  for  new  systems  because  they  are  often  empirical  in  nature)  or 
from  hydrocodes  (expensive  and  oftentimes  questioned) . In  Part  I a 
model  was  developed  to  predict  the  blast  field  which  was  simple  and 
economical  to  use  and  which  provided  insight  into  the  relevant  physi- 
cal processes  [1] . 

This  part  presents  an  overview  of  the  analysis  with  application 
to  guns  ranging  in  barrel  diameter  from  .30  caiber  to  8 inch. 


[1]  B.  Henriksen,  B.  Cummings, "A  Predictive  Scheme  for  the  Blast  En- 
vironment of  Army  Weapons,  Part  I.  Devlopment  and  Validation  of  the 
Theory,"  BRL  Technical  Report  ARBRL^TR-02044,  Feb.  1978.  (AD  #A053487) 
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II.  Theory 


The  approach  in  predicting  the  blast  overpressure  consists  of  two 
parts:  (1)  analysis  of  the  interior  losses  leading  to  a prediction  of 
an  equivalent  yield  for  a spherical  explosion  at  the  muzzle  and  (2) 
analysis  of  the  blast  wave  from  the  explosion  with  suitable  modifica- 
tions to  account  for  the  geometric  asymmetries  introduced  by  the  high 
exit  velocity  of  the  gases.  The  following  is  designed  to  give  an 
overview  of  the  analysis  — the  details  are  contained  in  Reference  1. 

The  analysis  is  an  accounting  of  the  energies  in  the  firing  pro- 
cess. Fundamental  to  the  theory  is  the  partitioning  of  energy  into 
either  prompt  energy  or  waste  heat.  Prompt  energy  is  that  energy 
which  is  "promptly"  available  for  supporting  the  shock  wave.  It  con- 
sists of  the  static  overpressure  and  the  dynamic  pressure.  Waste  heat 
includes  all  other  energies.  Examples  of  this  type  are:  energies 
which  for  some  reason  are  delayed;  energies  which  appear  as  ambient 
temperature  rise;  energies  which  perform  irreversible  wcrk  such  as 
crushing  a material, etc.  This  division  of  energies  is  important 
because  it  allows  focusing  the  analysis  only  on  the  prompt  energy  — 
once  an  energy  is  determined  to  be  waste  heat  it  can  be  excluded  from 
the  analysis. 

The  following  two  sections  address  the  interior  losses  and  the 
blast  wave  predictions  for  the  equivalent  explosion  at  the  muzzle. 


1 1 . 1 Interior  Analysis 

The  losses  from  the  initial  total  propellant  energy  are  grouped 
as  follows: 

a)  energy  imparted  to  the  projectile, 

b)  energy  lost  to  the  barrel  in  the  form  of  recoil,  barrel 
heating,  barrel  expansion,  etc., 

c)  energy  converted  to  waste  heat  by  formation  of  the  boun- 
dary layer,  and 

d)  energy  trapped  in  the  grooves  (rifling) . 

Energy  loss  a)  is  self-evident.  The  projectile  energy  does  not 
support  the  main  blast  wave.  A certain  amount  of  the  energy  does 
appear  as  a shock  wave  off  the  projectile;  however,  this  is  not  the 
strong,  damage-producing  wave. 

Energy  loss  b) : This  loss  is  estimated  from  experience  with  coni- 
cal shock  tubes.  Measurements  show  that  only  about  16%  of  the  initial 
explosive  energy  is  imparted  to  the  shock  wave.  This  can  be  visual- 
ized with  the  help  of  Figure  1.  If  we  view  the  propellant  as  a cube 
and  equipartition  the  direction  of  energy  flow  we  see  that  1/6  is 
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directed  towards  the  muzzle,  2/3  is  directed  to  the  walls  and  1/6  is 
directed  towards  the  rear  (recoil) . We  are  not  suggesting  that  the 
5/6  portion  of  the  energy  does  not  exit  the  tube,  only  that  it  is 
delayed  and  hence  is  not  promptly  available. 


RECOIL 

DIRECTION 


BARREL  AXIS 


Energy  loss  c) : We  propose  that  as  the  projectile  travels  down 
the  barrel  a boundary  layer  is  formed  as  depicted  in  Figure  2. 


The  boundary  layer  closure  divides  the  flow  into  two  distinct 
regions.  In  the  boundary  layer  the  flow  is  turbulent  with  the  energy 
partitioned  between  the  three  translational  degrees  of  freedom.  In 
the.  conical  region  between  the  closure  point  and  the  projectile  the 
motion  of  the  fluid  particles  is  ordered  and  traveling  in  the  direc- 
tion of  the  projectile.  We  propose  that  the  strength  of  the  blast 
wave  is  determined  by  the  pressure  in  the  fluid  at  the  projectile 
base.  In  the  turbulent  boundary  layer  a fluid  particle  has  a proba- 
bility of  traveling  in  any  of  the  three  directions;  hence  we  treat  this 
energy  as  delayed. 

We  can  estimate  the  pressure  at  the  projectile  base  by  treating 
the  ordered  flow  region  as  consisting  of  an  expansion  from  the  closure 
point  where  the  ordered  motion  is  essentially  zero  (in  a reference 
frame  associated  with  the  projectile)  to  some  velocity  at  the  projec- 
tile base.  TO  make  the  analysis  tractable  we  assume  there  are  no 
expansion  or  compression  waves  in  this  region  and  apply  the  approach 
of  Bernoulli  to  this  expansion  process. 
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Since  we  are  interested  in  motion  only  along  the  barrel  axis  (z 
axis)  the  appropriate  differential  equation  is 


(II. 1.1) 


0 


in  which  the  ordered  velocity  is  represented  by  u. 

In  Part  I an  examination  of  experimental  data  showed  that  the 
thermal  velocity  of  the  fluid  particles  in  random  motion  is  of  the 
same  order  as  the  projectile  velocity.  This  suggests  that  the  dis- 
tance from  the  closure  point  to  the  projectile  base  is  constant  with 
the  consequence  that  the  temporal  variation  of  u is  negligible. 

Equation  (II. 1.1)  becomes  an  ordinary  differential  equation  and, 
with  the  use  of  the  adiabatic  expansion  relations,  can  be  integrated 
to  yield 


(II. 1.2) 


constant 


where  Y is  the  ratio  of  the  specific  heats  and  a is  the  speed  of 
sound.  The  constant  is  evaluated  by  noting  that  at  the  closure  point 
the  ordered  velocity  is  essentially  zero  (stagnated)  and  hence  we 
associate  the  pressure  at  that  point  with  the  chamber  pressure. 

After  some  algebra  we  find  that  the  pressure  ratio  at  the  projec- 
tile base  can  be  determined  from 


(Pr  - l)2 


c 


(II. 1.3) 


P 


” <pr  +£r> 

P P 


r 


c 
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where  denotes  the  pressure  ratio  relative  to  ambient  and  the  sub- 
scripts c and  p refer  to  the  chamber  and  projectile  respectively.  For 
large  chamber  pressure  ratios  (>>1)  the  above  relation  becomes, 
approx imately , 


P = 18.9  P 

c p 


This  states  that  the  effect  of  the  choke  formation  is  a reduction  in 
the  prompt  energy  by  a factor  of  approximately  19. 

Knowledge  of  the  point  at  which  boundary  layer  closure  occurs  is 
of  importance  in  determining  subseqent  losses  (item  d) . Experimental 
results  using  shock  tubes  show  that  the  distance  down  the  barrel  at 
which  closure  occurs  can  be  estimated  from  the  empirical  relation  [2] 


L = 15 

D ~ H0’1 


where  H is  the  roughness  factor  and  equal  to  the  ratio  of  the  rough- 
ness height  (rifling)  to  the  unimpeded  diameter  of  the  tube. 

The  last  energy  loss  mechanism,  d,  is  determined  after  boundary 
layer  closure  has  occurred  and  continues  along  the  remainder  of  the 
barrel.  The  specific  process  is  visualized  using  Figure  3 in  which 
the  rifling  is  depicted  as  vertical  structures  approximating  baffles. 


Figure  3 

The  loss  mechanism  is  simply  an  impedance  of  the  energy  flow  due  to 
entrapment  within  the  rifling. 


[2]  F.B.  Porzel,  "Study  of  Shock  Impedance  Effects  in  a Rough  Walled 
Tunnel,"  IDA  Log  No.  HQ  68-7324,  Mar.,  1969  (AD  684790). 
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The  analysis  determines  the  change  in  the  total  energy.  It  is 
assumed  proportional  to  the  kinetic  energy  (per  unit  volume) , e , 
and  the  volume  subtended  by  the  average  roughness  of  the  barrel,  which 
itself  is  a product  of  the  rifling  height,  h,  perimeter  s and  the 
distance ,dL.  That  is. 


(II. 1.4) 


— a S h dL 


Employing  the  standard  relations  of  compressible  flow  — adia- 
batic, Rankine-Hugoniot,  etc.  — Equation  (II. 1.4)  can  be  integrated 
to  yield 


T1  ln(APr)  -i^ln(APr  + B) 


r 


constant  - 2 a ^ 


(II. 1.5) 


where 


and 


AP 


r 


is  the  overpressure  ratio 
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The  constant  can  be  eliminated  by  evaluating  the  expression 
between  two  points  along  the  barrel.  Specifically,  we  select  the 
overpressure  ratio  at  the  point  where  the  boundary  layer  closure 
occurs.  This  allows  us  to  determine  the  overpressure  ratio  at  the  muz- 
zle. Denoting  the  left-hand  side  of  Equation  (II. 1.5)  by  I,  we  obtain 
as  the  impedance  loss  between  the  two  points 


I 


1 


2 ct 


h h h 

D D 


(II. 1.6) 


Experimental  measurements  with  shock  tubes  show  that,  for  barrel 
lengths  of  interest,  the  proportionality  constant  is  equal  to 
unity. [2] 

In  the  present  analysis  the  kinetic  energy  of  the  moving  gas  at 
the  instant  the  projectile  leaves  the  barrel  is  included  in  the  hydro- 
dynamic  yield.  The  initial  prompt  energy  of  the  equivalent  explosion 
is 


Yo  = F [ (Mn e “ KE  ) 6 + i M V2  ] 

u P 2 p m 1 


vdiere  Mp  is  the  propellant  mass,  e is  the  specific  energy  of  the  pro- 
pellent, KE  is  the  projectile  kinetic  energy,  6 is  the  barrel  loss 
factor  as  determined  from  the  previous  discussion  and  Vm  is  the  muzzle 
velocity.  An  initial  value  of  the  explosion  radius  is  required  in 
addition  to  the  initial  prompt  energy  in  order  to  use  the  QZQ 
hypothesis  (discussed  in  the  next  section)  to  predict  the  blast  over- 
pressures. For  the  gun  muzzle  blast,  the  initial  explosion  may  be 
considered  as  the  mass  of  high  pressure  gas  which  exits  the  muzzle 
immediately  after  the  projectile  has  left  the  muzzle.  Schlieren  pho- 
tographs taken  at  these  times  show  a "barrel  shock"  system  whose 
characteristic  size  is  of  the  order  of  two  muzzle  diameters.  Hence,  in 
the  present  analysis  the  barrel  diameter  is  taken  as  the  initial 
radius.  Further,  the  QZQ  scaling  is  relatively  insensitive  to  the 
initial  radius  chosen. 
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II. 2 Blast  Wave  Analysis 


The  core  of  the  analysis  for  the  blast  field  is  the  Unified 
Theory  of  Explosions  as  developed  by  F.B.Porzel . [3]  In  essence,  the 
analysis  is  simply  an  accounting  of  the  production  of  waste  heat  — 
the  remainder  of  the  initial  energy  being  in  the  blast  wave.  The  com- 
putational procedure  can  be  illustrated  using  Figure  4.  Figure  4 
shows  the  standard  Pv  diagram  for  an  expansion  process. 


The  compressed  gas  starts  at  some  initial  state,  Pv.  In  expand- 
ing back  to  ambient  the  gas  follows  the  curved  line,  the  so-called 
Hugoniot  or  shock  adiabat.  The  straight  line  connecting  the  initial 
and  final  states  represents  the  thermodynamic  path  a wave  of  unchang- 
ing shape  would  follow.  The  trajectory  below  the  straight  line  simply 
states  that  compression  waves  steepen  while  rarefaction  waves  expand. 

From  the  Rankine-Hugoniot  relations,  we  find  that  the  rectangle 
bounded  by  P = constant,  v = constant  and  P = 0,  v = constant  is  the 
total  energy:  a 


eT  = P (va  - v)  . 


(II. 2.1) 


The  upper  triangle  in  the  rectangle  is  the  kinetic  energy: 


eKE  = \ (P  " pa}  (va“v)*  (II. 2. 2) 


a Unified  Theory  of  Explosions 


[3]  F.B.  Porzel,  "Introduction  to 
(UTE),"  NOLTR  72-209,  Sept. ,1972 
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The  remaining  trapezoidal  area  is  the  internal  energy: 


e 


(II. 2. 3) 


The  area  between  the  adiabat  and  the  straight  line  connecting  Pv  and 
ambient  is  the  waste  heat.  At  an  infinite  distance  from  the  source  of 
the  blast  all  of  the  prompt  energy  has  been  converted  to  waste  heat. 

We  can  estimate  the  ratio  of  the  production  of  waste  heat  to  the 
total  energy  using  the  Rankine-Hugoniot  relations  for  changes  across  a 
shock  wave: 


1 <P  + Pa>  <va  - v> 


P + P 

a 


p (va  - V) 


2 P * 


For  P»P  the  internal  energy  is  approximately  1/2  the  total;  hence, 
half  the  energy  is  subject  to  waste  (the  kinetic  energy  is  not  subject 
to  waste) . In  the  acoustic  approximation,  P ~ P and  all  the  energy  is 
wasted . a 

In  the  following,  two  problems  are  addressed: 

a)  properties  of  the  spherical  blast  wave,  and 

b)  modifications  to  include  the  effect  of  the  high  direc- 
tionality of  the  exiting  gases. 


II. 2. a The  Unified  Theory  of  Explosions 


The  shock  expansion  process  is  adiabatic;  hence,  in  general  we  can 

write 


e. 


T 


W + eKE  + Q 


(II. 2. 4) 
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is  the  kinetic 


where  W is  the  pressure  volume  energy,  W = / Pdv,  e^ 
energy  per  unit  mass  and  Q is  the  specific  waste  heat. 

The  total  prompt  energy,  the  integral  of  the  prompt  energy,  is 
defined  by 


R 

(W  + e^,)  r2  dr.  (II. 2. 5) 

0 


Subject  to  the  boundary  conditions: 

a)  at  R = Ro,  the  initial  charge  radius,  Y is  the  hydro- 
dynamic  yield,  and 

b)  as  R approaches  infinity  the  shock  wave  must  completely 
dissipate, 

the  expression  for  the  total  prompt  energy  becomes 


Y(R) 


dr. 


(II. 2. 6) 


The  abstraction  which  makes  UTE  tractable  is  the  QZQ  hypothesis 
which  states: 


Q Z^  = constant. 


(II. 2. 7) 


where  Z is  a mass  corrected  radius  and  q is  a constant  (which  has 
values  of  3.5  in  the  strong  shock  regime  and  4.0  in  the  weak) . 
Specification  of  a relation  between  R and  the  mass  corrected  radius 
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permits  a closed  form  integration  of  Equation  (II. 2. 6).  The  proof  of 
Equation  (II. 2. 7)  is  involved  and  will  not  be  presented  here.  The 
reader  is  referred  to  either  [1]  or  [3] . 


The  concept  of  a mass  correction  to  the  radius  was  originally 
introduced  to  permit  analysis  of  entirely  different  blast  wave  sources 
— the  "massless"  nuclear  and  the  standard  TNT  charge.  The  relation 
between  the  two  radii  is  simply 


Z 


(II. 2. 8) 


where  M represents  (essentially)  the  ratio  of  the  energy  in  the  mass 
to  the  energy  in  the  air  contained  within  the  shock  boundary. 

The  only  remaining  piece  of  information  necessary  for  solution  of 
the  blast  wave  properties  is  determination  of  the  constant  in  Equation 
(II. 2. 7).  This  is  found  by  specifying  the  pressure  at  the  transition 
point  between  the  strong  and  weak  shock  regimes.  In  this  analysis  we 
have  selected  a pressure  ratio  of  2 since  this  is  approximately  the 
point  at  which  the  negative  phase  develops. [3]  Using  this  value  we  can 
solve  Equation  (II. 2. 7)  and  (II. 2. 8)  to  determine  the  mass  corrected 
transition  radius  and  finally  determine  the  constant  via 


Q Zg  = constant  = Qt  Z^. 


The  computational  loop  is  closed  by  determining  the  area  between  the 
Hugoniot  and  the  straight  line  (the  waste  heat  in  Figure  4).  It  is 
found  to  be 


Q 


* 


(II. 2. 9) 
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where  the  density  ratio  is  obtained  from  the  Rankine-Hugoniot  rela- 
tions : 


D 


P 


P 


a 


(II. 2. 10) 


II. 2. b Asymmetric  Blast  Effects 


Experimental  measurements  of  the  blast  fields  produced  by  guns 
have  shown  that  for  a fixed  distance  from  the  observer  to  the  muzzle, 
considerably  larger  overpressures  occur  in  the  region  forward  of  the 
muzzle  than  in  the  region  behind  the  muzzle.  This  is  analogous  to  the 
moving  charge  effect  which  has  been  studied  by  Armendt  and  Sperrazza 
[4] . They  found  that  the  blast  from  a moving  charge  remained  essen- 
tially spherical  about  an  origin  which  moved  with  the  center  of  mass 
of  the  decelerating  charge.  The  deceleration  of  the  charge  could  be 
determined  by  conservation  of  momentum.  The  center  of  mass  decelerates 
rapidly  as  the  expanding  shock  wave  engulfs  an  ever-growing  mass  of 
initially  stationary  air. 

Detailed  calculations  were  made  based  on  this  effect  and  while  it 
predicts  higher  overpressures  forward  of  the  muzzle  than  aft,  the 
effect  is  too  small  to  explain  the  order  of  magnitude  differences 
which  are  observed  in  the  muzzle  blast  experiments.  This  is  due,  pri- 
marily, to  the  rapid  deceleration  of  the  center  of  mass  which  is  pro- 
duced by  a rapidly  expanding  spherical  shock. 

Porzel  [3]  described  a concept  called  generalized  divergence 
(GDV)  which  permits  an  extension  of  UTE  to  non-spherical  geometries. 
He  notes  that  the  physical  significance  of  the  spatial  coordinate  in 
the  hydrodynamic  equations  is  a radius  of  curvature  rather  than  the 
location  relative  to  seme  earlier  position.  It  is  the  local  radius  of 
curvature  of  the  wave  front  which  determines  its  divergence.  Thus  it 
is  the  initial  shape  of  the  charge  which  determines  the  shape  of  the 
blast  field. 

[4]  B.F.  Armendt,  J.  Sperrazza,  "Air  Blast  Measurements  Around  Moving 
Explosive  Charges,  Part  III,"  BRL  Memorandum  Report  No.  1019,  July, 
1956  (AD  #114950) 
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The  initial  charge  for  a gun  muzzle  blast  is  the  barrel  gases 
which  exit  when  the  projectile  leaves  the  muzzle.  The  high  overpres- 
sure of  these  gases  causes  a "barrel  shock"  system  to  form  as  shown 
schematically  in  Figure  5.  This  provides  an  initial  source  for  the 
blast  field  which  is  far  from  spherical.  The  local  radius  of  curva- 
ture is  much  greater  in  the  forward  portion  of  the  shock  system, 
resulting  in  a smaller  divergence  and  less  rapid  decrease  in  the  blast 
overpressure  forward  of  the  gun  muzzle  compared  to  the  rear  where  the 
radius  of  curvature  is  much  smaller. 

The  concept  of  GDV  was  incorporated  into  the  scaling  of  blast 
overpressure  provided  by  the  QZQ  hypothesis  by  modifying  the  distance 
scale  at  each  angular  position  in  the  blast  field  to  correct  for  the 
initial  non-spherical  geometry.  The  value  of  R used  in  the  QZQ  calcu- 
lation was  given  by 


(II. 2.11) 


where  G(e)  is  a geometry  factor  which  was  determined  on  an  ad  hoc 
basis.  It  was  found  that  the  experimental  data  were  reasonably  well 
represented  using  a geonetry  factor  given  by 


(II. 2. 12) 


The  derivation  of  this  function  was  based  on  a velocity  argument  and 
is  given  in  Appendix  A. 


II. 3 Pulse  Length  Calculations 


In  Part  1 the  pulse  length  of  the  overpressure  pulse  was  deter- 
mined based  on  a characteristic  length  and  the  particle  velocity.  The 
characteristic  length  was  representative  of  the  volume  occupied  by  the 
remaining  prompt  energy  in  the  blast. 

An  alternative  method  of  calculating  the  pulse  length  has  been 
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FIGURE  5.  "Barrel  Shock"  Produced  by  Propellent  Gases  Exiting  Muzzle. 
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developed  based  on  the  propagation  of  waves  of  finite  amplitude.  As 
the  blast  wave  from  the  explosion  expands,  the  initial  high  overpres- 
sure at  the  shock  front  decreases  behind  the  shock  front.  At  some 
point  during  the  expansion  the  overpressure  behind  the  shock  drops  to 
zero  and  upon  further  expansion  the  blast  wave  develops  a negative 
phase. [5]  The  radius  at  which  the  overpressure  first  drops  to  zero  is 
called  the  transitions  radius.  The  fluid  velocity  behind  the  shock 
has  a similar  behavior  with  a large  velocity  in  the  direction  of  the 
shock  velocity  decreasing  to  zero  and  becoming  negative.  The  result- 
ing pulse  forms  are  shown  in  Figure  6.  These  pulses  of  finite  ampli- 
tude propagate  into  the  undisturbed  air.  The  propagation  velocity  is 
not  the  same  for  all  portions  of  the  pulse,  but  depends  on  the  local 
speed  of  sound  and  the  fluid  velocity.  The  leading  edge  of  the  pulse 
propagates  with  a greater  velocity  than  the  rest  of  the  pulse  and 
therefore  the  pulse  increases  in  length  as  it  propagates  outward. 

The  local  wave  speed  for  the  pulse  of  finite  amplitude  is  given 
by  [6] 


c = a + u 
n 


(II. 3.1) 


where  a and  u are  the  local  value  of  the  speed  of  sound  and  fluid 
velocity.  For  an  isentropic  process  the  local  wave  speed  becomes 


c = 


(II. 3. 2) 


where  a is  the  speed  of  sound  in  the  undisturbed  region  ahead  of  the 
pulse.  The  propagation  of  a shock  is  not  an  isentropic  process,  but 
for  overpressures  at  distances  greater  than  the  transition  radius  the 
corrections  are  small  compared  to  the  uncertainty  in  the  measured 
values  of  pulse  length. 

The  trajectories  of  the  pulse  front  and  the  point  at  which  the 
fluid  velocity  drops  to  zero  are  shown  in  Figure  7.  The  pulse  shape 
is  shown  at  the  transition  radius  R^_  and  at  two  other  positions. 

We  will  define  the  pulse  length  t as  the  time  between  these  two 
trajectories  at  a given  position  R.  Note  that  as  the  shock  wave 


[5]  Yu.B.  Zel'dovich,  Yu.P.  Raizer,  Physics  of  Shock  Waves  and  High 
Temperature  Hydrodynamic  Phenomena,  Vol.  1,  Academic  Press,  New  York, 
1966. 

[6]  H.W.  Liepmann,  A.  Roshko,  Elements  of  Gasdynamics,  John  Wiley  & 
Sons,  New  York,  1957. 
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Velocity  Overpressure 
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FIGURE  6. 
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Sketch  of  Overpressure  and  Fluid  Velocity  in  the  Blast  Pulse 
After  the  Negative  Phase  Has  Developed. 
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* 


FIGURE  7.  Sketch  of  the  Trajectories  of  the  Pulse  Front  and  the 
Point  Where  Fluid  Velocity  Drops  to  Zero. 
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expands  and  the  fluid  velocity  behind  the  shock  approaches  zero,  the 
shock  front  approaches  the  speed  of  sound  and  the  pulse  length 
approaches  an  asymptotic  value  which  then  propagates  as  an  acoustic 
wave. 


The  zero-velocity  trajectory  is  given  by 


+ + xt 


(II. 3. 3) 


and  the  shock  front  velocity  by 


t 


2 


dx 


1 + 


TT * + t.  . 

+1  U t 


(II. 3.4) 


The  difference  of  these  expressions  gives  the  pulse  length 


T 


Tt  + i 


<R-Rt> 


dx 


+ I±i  H 

2 a 


(II. 3. 5) 


where  is  the  initial  pulse  length  at  the  transition  radius  R.  . The 
initial  pulse  length  can  be  estimated  with  the  aid  of  Figure  8 
which  shows  the  overpressure  at  the  time  t when  the  shock  reaches  the 
transition  radius. [7]  This  is  the  radius  for  which  the  negative  phase 
has  fully  developed  and  occurs  at  a pressure  ratio  of  two  across  the 

m H.L.  Brode,  "Numerical  Solutions  of  Spherical  Blast  Waves,"  J. 
Appl.  Phys.,  Vol.  26,  No.  6,  June,  1955,  pp.  766-775 
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FIGURE  8. 


Overpressure  Pulse  at  Transition.  The  Pressure  Ratio 
Drops  to  Unity  of  Approximately  One-Half  the  Transition 
Radius  r^ . 
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shock.  If  we  make  the  ad  hoc  assumption  that  the  pulse  occupies  half 
the  spherical  radius  at  transition  then  the  characteristic  initial 
pulse  length  is  given  by 


(II. 3. 6) 


where  c is  the  wave  speed  of  the  shock  at  a pressure  ratio  of  two. 
s 

The  pulse  length  is  calculated  from  Equation  (II. 3. 5)  using  the 
initial  pulse  length  given  in  Equation  (II. 3. 6)  and  the  fluid  velocity 
determined  frcxn  the  Rankine— Hugoniot  relations  and  the  local  value  of 
overpressure  determined  frcm  UTE. 


III.  Results 

The  purpose  of  this  study  was  to  determine  the  applicability  of 
the  theory  presented  in  Part  I for  smaller  caliber  guns  and  to  extend 
the  analysis  to  include  those  positions  aft  of  the  muzzle  exit  plane. 


III.l  Guns  Considered 

Ihe  above  theories  were  applied  to  six  different  guns  and  to  one 
case  where  the  propellent  mass  and  specific  energy  were  changed  (zone 
5 vs.  zone  8) . The  guns  and  their  characteristics  are  listed  in  Table 
1.  [8] [9]  The  barrel  diameters  range  frcm  the  .30  caliber  pistol/rifle 
to  an  8 inch  naval  gun. 

The  .30  caliber  pistol  was  (apparently)  a shortened  .30  caliber 
rifle.  All  of  the  physical  characteristics  are  the  same,  the  different 
projectile  muzzle  velocity  being  obtained  from  a different  propellant 
loading.  The  two  cases  using  the  105  mm  Howitzer  provide  an  excellent 
evaluation  of  the  effect  of  muzzle  velocity  on  the  initial  yield. 


[8J  P.S."  "Wes  tine,  J.C.  HoKanson,  "Prediction  of  Stand-off  Distances  to 
Prevent  Loss  of  Hearing  frcm  Muzzle  Blast,"  Southwest  Res.  Inst.  Re- 
port No.  R-CR-75-003,  Feb.,  1975  (ADA005274). 

[9]  P S.  Westine,  "Modeling  the  Blast  Field  Around  Naval  Guns  and  Con 
ceptual  Design  of  a Model  Gun  Blast  Facility,"  Southwest  Res.  Inst. 
Report  No.  TR  02-2643-01,  Sept.,  1970  (AD  875984). 
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All  the  chamber  pressure  ratios  were  assumed  approximately  3000  — the  results  are  very  insensitive 
to  this  parameter. 


III. 2 Comparison  of  Theory  and  Experiment 


The  predictions  of  the  theory  have  been  compared  to  these  data 
and  the  results  are  presented  in  Figures  9 through  13.  The  theory  is 
clearly  capable  of  predicting  the  muzzle  blast  overpressure  for  guns 
of  vastly  different  calibers.  Careful  examination  of  the  curves  shows 
that  the  theory  will  also  predict  the  angular  variations  in  the  blast 
field  with  reasonable  accuracy  over  the  full  range  of  calibers  exam- 
ined. The  largest  discrepency  between  theory  and  experiment  occurs  at 
distances  beyond  50  calibers  where  the  experimental  data  exceeds  the 
theoretical  predictions.  The  experimental  data  were  all  obtained  with 
the  guns  firing  essentially  horizontally  over  the  ground  or  a ground 
plane.  The  shock  wave  associated  with  the  muzzle  blast  will  reflect 
from  this  plane  causing  a reinforcement  to  the  primary  shock  which 
increases  the  measured  overpressure.  Indeed,  in  much  of  the  data  the 
measured  overpressure  values  increase  with  increasing  distance  from 
the  muzzle  for  distances  beyond  50  calibers.  Since  the  theory  does  not 
account  for  reflected  shocks  the  disagreement  at  large  distances  is 
not  surprising. 

Examination  of  the  angular  behavior  of  the  theory  for  distances 
less  than  50  calibers  for  the  cases  considered  reveals  that  on  the 
average  the  theoretical  predictions  tend  to  overestimate  the  overpres- 
sure in  the  aft  (greater  than  90  degrees)  portion  of  the  blast  field 
and  underestimate  the  overpressure  in  the  forward  portion  of  t"he 
field.  This  is  particularly  noticeable  for  the  .30  caliber  pistol, 
the  smallest  gun  considered.  This  effect  is  a result  of  the  particu- 
lar form  of  the  geometry  factor  chosen.  An  improvement  might  be 
obtained  by  modifying  the  form  of  the  geometry  factor  according  to  the 
ratio  of  random  to  ordered  motion  in  the  muzzle  gases  as  proposed  in 
Part  I. 

Since  the  current  theory  contains  the  kinetic  energy  of  the  muz- 
zle gas  in  the  prompt  energy  the  overpressure  was  calculated  for  a 105 
mm  howitzer  in  which  the  muzzle  velocity  varied  as  a result  of  dif- 
ferent propellent  charges.  The  results  are  shown  in  Figure  14.  The 
theory  accurately  predicts  muzzle  blast  overpressure  for  these  cases 
although  the  muzzle  velocity  is  higher  by  a factor  of  two  for  the  zone 
8 (Z8)  compared  to  the  zone  5 (Z5)  experiment.  The  result  shows  that 
the  theory  will  predict  the  muzzle  blast  from  guns  over  a considerable 
range  of  muzzle  velocity. 

Pulse  length  of  the  overpressure  pulse  is  more  difficult  to  com- 
pare since  there  are  large  variations  in  the  experimental  values 
obtained.  A typical  situation  is  shown  in  Figure  15  for  the  case  of 
the  same  20  mm  gun  shown  in  Figure  12.  The  theory  predicts  that  the 
pulse  length  will  increase  with  distance  from  the  muzzle  until  it 
reaches  some  asymptotic  value  but  the  experiments  often  show  that  the 
pulse  length  first  increases  then  decreases  with  distance.  This 
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Radial  Distance  From  Muzzle  - Calibers 

FIGURE  15.  Overpressure  Pulse  Length  from  20  mm  Gun. 
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Radial  Distance  From  Muzzle  - Calibers 
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Radial  Distance  From  Muzzle  - Calibers 
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Radial  Distance  From  Muzzle  - Calibers 
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FIGURE  14.  Overpressure  from  105  mm  Howitzer  at  Two  Muzzle 
Velocities.  ZS  is  650  m/s  and  Z5  is  332  m/ s . 


decrease  first  appears  at  distances  greater  than  50  calibers  and  may 
be  related  to  the  shock  reflection  from  the  ground  plane  which  pro- 
duces pulses  of  complex  shape. 

Predicted  pulse  lengths  for  four  guns  ranging  from  an  8 inch 
naval  gun  to  a .30  caliber  rifle  are  shown  in  Figure  16.  The  pulse 
length  was  calculated  for  an  angle  of  90  degees  to  the  line  of  fire 
and  the  experimental  data  were  taken  at  75  and  105  degrees.  The  meas- 
ured values  of  the  pulse  length  are  reasonably  well  predicted  by  the 
theory  over  the  wide  range  of  pulse  lengths  produced  by  guns  of 
greatly  different  caliber. 

These  results  indicate  that  the  theory  gives  reasonable  agreement 
with  experiment  over  a wide  range  of  gun  calibers.  The  theory  has  the 
advantage  of  simplicity  and  ease  of  calculation  and  should  be  of  great 
value  in  the  prediction  of  free-field  gun  muzzle  blast  effects. 
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FIGURE  16.  Overpressure  Pulse  Length  for  Four  Guns. 
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APPENDIX 


The  form  of  the  geometry  function  G(9)  was  developed  on  the  basis 
of  a moving  charge.  If  the  velocity  at  some  point  in  the  blast  field 
is  the  sum  of  the  charge  veloctiy  V and  the  blast  velocity  for  a sta- 


tionary charge  V then  V = V + Vc 
of  the  radius  vector  from  the  Suzzli 


is  assumed  to  be  in  the  direction 
If  the  field  position  is  deter- 
mined by  the  distance  from  the  muzzle  and  the  polar  angle  relative  to 
the  line  of  fire  then 


0 is  the  angle  between  V and  V and 

c 


V + V 


or 


i-Y  - 2 


VsVs 


V2 

c 


cos  e + 


- 2WcCOS  0 


V 


c\2 


V. 


-1  = 0. 


Solving  for  the  ratio  V/Vs  we  obtain 


_V 

V. 


V 
c 

V 


COS  9 


t J/£ 


/V  \„" 

/ 2 

i / s\ 2 

cos  e - 

1 - — 

L ' C'  -i 

s s 

We  assume  that  the  geometry  factor  has  a similar  form,  i.e.. 


G(9)  = a | cose 

where  <*  and  0 are  constants  to  be  determined 


cos2e 


+ 6‘ 


(Al) 


(A2) 


(A3) 


(A4) 


If  we  require  that  at  an  angle  of  90  degrees  from  the  barrel  axis 
the  geometry  factor  have  the  value  unity,  then 


3 = 1/n 


In  Part  1 [1]  we  found  that  the  kinetic  energy  of  ordered  motion  at 
the  muzzle  was  approximately  twice  that  of  random  motion.  Based  on 
this  result  we  require 


G(0)  = 2. 


This  gives  the  results 


a = 


3 

4 


and 


G(0) 


3 

¥ 


cos  e + 


V' 


cos2e 


(A6) 


(A7) 


(A8) 
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BASIC  PROGRAM 


A BASIC  computer  program  has  been  constructed  to  perform  the  cal- 
culations for  the  muzzle  blast  field  based  on  the  preceding  analysis. 
The  program  is  listed  at  the  end  of  this  section.  The  input  parameters 
required  for  the  program  are  as  follows: 

1.  Projectile  mass  - kilograms  (kg) 

2.  Muzzle  velocity  - metres/second  (m/s) 

3.  Propellant  specific  energy  - joules/kilogram  (j/kg) 

4.  Ratio  of  peak  chamber  pressure  to  ambient  - dimensionless 

5.  Barrel  length  - metres  (m) 

6.  Barrel  diameter  - metres  (m) 

7.  Height  of  barrel  rifling  - metres  (m) 


When  execution  of  the  program  begins  the  program  will  pause  for 
input  until  these  values  are  entered  in  the  order  given  above.  When 
the  last  value  is  entered  the  program  calculates  the  overpressure  at 
the  choke  position  and  at  the  muzzle  and  an  equivalent  explosive 
yield.  Next,  the  program  will  pause  to  allow  optional  values  of  Ql, 
Q2  and  the  ratio  of  specific  heats  to  be  entered  if  desired.  The 
variables  Ql  and  Q2  are  the  exponents  in  the  Porzel  QZQ  hypothesis  for 
the  strong  and  weak  shock  regimes  and  have  the  standard  values  of 
3.25  and  3.5  respectively.  The  standard  value  of  gamma  is  1.4.  After 
the  option  is  exercised  the  program  computes  the  transition  radius  and 
pauses  until  the  desired  angular  position  in  the  blast  field  is 
entered.  This  is  the  angle  measured  from  the  line  of  fire  to  the 
radius  vector  from  the  muzzle  to  the  field  location  and  lies  between  0 
and  180  degrees. 

Hie  program  is  designed  to  calculate  overpressure  and  pulse 
length  as  a function  of  distance  from  the  muzzle  at  the  input  angle. 
The  program  will  pause  until  the  desired  distance  closest  to  the  muz- 
zle is  entered.  It  then  pauses  until  an  increment  in  the  radial  dis- 
tance is  entered  and  again  until  the  total  number  of  increments 
desired  is  entered.  The  program  then  computes  the  overpressure  and 
pulse  lengths  at  each  radial  location  specified  for  the  input  angle. 
When  these  computations  are  completed  the  program  pauses  to  allow  the 
option  of  terminating  the  program  or  specifying  another  angle.  If  a 
new  angle  is  entered  then  the  program  pauses  to  allow  an  option  of 
continuing  with  the  same  radial  locations  or  changing  to  new  radial 
locations.  When  the  last  angular  position  desired  has  been  calculated 
the  operator  enters  minus  one  to  terminate  execution. 

Please  note  that  this  listing  is  in  lower  case  and  a mix  between 
BASIC  and  BASIC+.  The  user  must  change  it  for  his  particular  systan. 
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PROGRAM  LISTING 


100 

110 

120 

140 

150 

for  n=l  to  5 step  1 

print 

next  n 

print  "Muzzle  Blast  Overpressure  and  Pulse  Length' 
pr  in  t 11 

160 

print 

170 

print 

180 

print 

190 

rem. . . input  section 

200 

rem. . .variable  definitions  are  as  follows 

201 

rem. . .projectile: 

202 

rem, . . 

203 

rem. . .propellant: 

204 

rem, . . 

205 

rem. . . 

206 

rem. . .pressure: 

207 

rem. . .barrel 

208 

rem. . . 

209 

rem. . . 

210 

print  "projectile" /’mass  (kg)",, 

220 

input  wO 

225 

print 

230 

print  /’velocity  (m/sec)". 

240 

input  vO 

245 

print 

246 

print 

250 

print  "propellant" , "mass  (kg)",. 

260 

input  wl 

265 

print 

290 

print, "specific  energy  (j/kg)". 

300 

input  x2 

305 

print 

306 

print 

310 

print  "chamber  pressure  ratio  (p/p0)",. 

320 

input  p0 

325 

print 

326 

print 

330 

print  "barrel", "length  (m)",. 

340 

input  L0 

345 

print 

350 

print  , "d iameter  (m)  " , , 

360 

input  dO 

365 

print 

370 

print  /’groove  height  (m) " 
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380  input  LI 
390  print 
395  print 

400  rem. . . end  input  data 
1000  rem... 

1010  rem... begin  energy  calculations 

1011  rem. . . 

1012  rem... 

1013  rem. . . 

1020  e2  = x2*wl 
1030  el  = .5*w0*v0~2 
1040  eO  = e2  - el 
1050  print 

1060  print  energy  calculations" 

1070  print  n 

1080  print  "total  energy  in  propellant" , ,e2; "joules 
1090  print  "kinetic  energy  in  projectile" , ,el; "joules" 

1100  print  "maximum  available  energy" , ,e0; "joules" 

1110  pO  = p0  - 1 

2000  print 

2001  print  interior  losses" 

2002  print 

2010  rem...  roughness  factor  = h,  overpressure  at  choke  - x 

2020  rem...  o'pressure  ratio  at  muzzle  and  initial  yield  = y0 

2021  rem...  choke  length  = L2 
2030  h = Ll/dO 

2040  L2  = 0 

2050  x = p0  , , 

2060  rem...  if  choke  1/d  > gun  then  1/d  assume  no  choke  losses 

2065  if  h = 0 goto  2180 

2070  if  15/h~ . 1 > (LO/dO)  goto  2180 

2080  L2  = 15/h~.l 

2090  rem...  if  p0  > 100  use  fitted  curve 
2100  if  p0>100  goto  2160 
2110  x = p0 

2120  c = (l+x)*(l+4*x~2/5/(l+x)/(9+x))  5-1 
2130  if  abs(c-p0)<.001  goto  2170 
2140  x = x - (x-1) * (c-p0)/ (c-1.29) 

2150  goto  2110 

2160  x = exp( . 9211*log (pO)  - 2.138) 

2170  print  "overpressure  ratio  at";L2*dO;"m  is" ,x;  (choke) 
2180  y0  = exp( (2*h* (L2-L0/d0)+1.068*log (x) ) /1.068) 

2190  print  "overpressure  at";L0;"  m is",,y0;"  muzzle" 

2200  f = p0/y0 

2210  rem...  reduction  in  e (avail) 

2220  y0  = e0/(6*f) 

2230  e3  = wl/12*v0~2 

2231  yO  = y0  + e3 

2232  print  . ^ ,, 

2233  print  ,"  equivalent  explosion  parameters 
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2234 

2240 

2250 

2251 

2252 

2253 

2254 

2260 

2261 

2262 

2270 

2271 

2272 

2273 

2274 

2275 

2276 

2280 

2281 

2282 

2300 

2301 

2302 

2303 

2304 

2305 

2310 

2320 

2400 

2420 

2430 

2440 

2450 

2451 

2452 

2453 

2454 

2455 

2456 

2460 

2461 

2480 

2481 

2482 

2500 

2510 

2520 

2530 

2540 

2600 


print 

print  "yield",,, y0;"  joules" 
rO  = dO 
mO  = wl/6 

for  m = 1 to  5 step  1 

print 

next  m 

print  "for  optional  ql,  q2,  gamma  enter  one;  std.  values  enter  0" 
input  ol 

if  ol  = 0 goto  2280 
print  "enter  ql" 
input  ql 

print  "enter  q2" 
input  q2 

print  "enter  ratio  of  specific  heats,  gamma" 

input  g 

goto  2300 

ql  = 3.25 

q2  = 3.5 

g = 1.4 

rem...  qO  = pressure/density  (atmospheric) 

q0  = 78340 

h = .25 

q5  = 1889 

p0  = 101325 

d = 1.293 

m = m0*h/5.416 

z0  = (r0~3  + mr(l/3) 

rem...  find  transition  radius 

al  = (ql-q2)/(3-q2) 

a2  = (ql-3)*y0/(12.566*q5*z0~3) 

if  ql=q2  goto  2480 

z = 1 

s = z 

f = al*s~3  - s~ql  + a2 
fl  = 3*al*s~2  - ql*s~(ql-l) 
z = z - f/fl 

if  abs  (z-s)  < .001  goto  2460 

goto  2451 

zl  = z 

goto  2481 

zl  = a2~ (1/ql) 

r2  = ( (zl*z0)  ~3  - mr(l/3) 

print  "transition  radius" ,r2, "metres" 

rem...  set  up  qzq  constants  a and  b 

a = q5/q0*(zl/'ql) 

b = q5/q0*  (zl/'q2) 

print 

print 

rem. . . set  up  desired  radial  locations 
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2601  print  "input  angle  (degrees)" 

2602  input  tl 

2603  let  t2=cos(tl*pi/180) 

2604  r8=.75*(t2  + sqr (t2"2+16/9) ) 

2606  print  "enter  initial  radius  (metres)" 

2607  input  r3 

2610  print  "input  radial  increment  (metres)" 

2611  input  r4 

2620  print  "input  number  of  increments  desired" 

2621  input  k 

2622  print  "radius" /’overpressure" /’pulse  length" 

2623  print  " (metres) " (p-p0)/p0" ( seconds) " 

2624  print  "********","***********","************" 

2625  si  = 0 

2626  cO  = sqr(g*p0/d) 

2627  t5  = r2*r8/1000 

2630  for  i=l  to  k 

2640  r = (r3+r4*(i-l) ) 

2650  for  j=l  to  10 

2660  if  r <=  r2*r8  goto  2780 

2670  if  (r-r2*r8)  < r4  goto  2700 

2680  y2  = r4/10 

2690  goto  2710 

2700  y2  = (r-r2*r8)/10 

2710  v = (r-(10-j) *y2)/r8 

2720  gosub  4000 

2730  u=  (pl-l)*sqr(2*q0/((g+l)*pl+(g-l))) 

2740  si  = si  + y2/(l+(g+l)/2*u/c0) 

2750  next  j 

2760  t = t5  + ( (r-r2*r8)  - sl)/c0 

2770  goto  2800 

2780  t = t5 

2785  v = r/r8 

2790  gosub  4000 

2800  p = pl-1 

2810  print  r,p,t 

2820  goto  2840 

2830  print  "error;  radius"  ,r /’less  than  muzzle  diameter’’ 
2840  next  i 

3040  print  "input  next  angle  (degrees)  or  -1  to  terminate 

3041  input  tl 

3050  if  tl<0  goto  3070 

3051  t2  = cos(tl*pi/180) 

3052  r8  = .75*(t2+sqr (t2"2+16/9) ) 

3060  print  "enter  1 to  change  radius;  0 to  continue" 

3061  input  a8 

3062  if  a8  = 0 goto  2625 

3063  goto  2606 

3070  print  "end  of  ute  muzzle  blast  calculations" 

3080  goto  9999 
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4000  if  v<=r0  goto  2830 

4010  z = ((v‘3-Hn)<‘(l/3))/zO 

4020  if  z <=  zl  goto  4060 

4030  q = b/(z''q2) 

4040  yl  = q0*q/(q2-3)*(z0*z)''3 
4050  goto  4070 

4070  rem. . . now  calculate  overpressure 
4080  if  q < 7.45e-7  goto  4210 
4090  if  q > 1.167  goto  4310 

4100  pi  = 5 

4110  al  = (g+l)/(g-l)  + pi 

4120  a2  = (g+l)/(g-l)*pl  + 1 

4130  a3  = l+(g-l)*q 

4140  f = al*pl"(l/g)  - a3*a2 

4150  fl  = l/g*al*pl/'  (1/g-l)  + pl"(l/g) 

4160  fl  = fl  - (g+l)/(g-l) *a3 

4170  n = f/fl 

4180  if  abs(n)  < .001  goto  4340 
4190  pi  = pl-n 

4200  goto  4110 

4210  p=.02 

4220  cl  = 24*q*g~3/(g+l) 

4230  f = 3*p~4  - 2*p"3  + cl 

4240  fl  = 12*p~3  - 6*p~2 

4250  n = f/fl 

4260  if  abs(n)  < .0005  goto  4290 

4270  p = p-n 

4280  goto  4230 

4290  pi  = ptl 

4300  goto  4340 

4310  kl  = 22+16*log (q) /2 . 303 

4320  1 = 11.5  - sqr (132.25  - kl) 

4330  pi  « 10*1  + 1 

4340  return 

9999  end 
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P.  0.  Box  482 
Fort  Worth,  TX  76101 
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Copies  Organization  Copies  Organization 


1 Bell  Telephone  Labs,  Inc. 

ATTN:  Tech  Rpt  Ctr 
Mountain  Avenue 
Murray  Hill,  NJ  07971 

1 The  Boeing  Company 
ATTN:  Aerospace  Library 
P.  0.  Box  3707 
Seattle,  WA  98124 

2 California  Research  and 

Technology,  Inc. 

ATTN : Ken  Kreyenhagen 

Technical  Library 
6269  Variel  Avenue 
Woodland  Hills,  CA  91364 

1 Calspan  Corporation 

ATTN:  Technical  Library 
P.  0.  Box  235 
Buffalo,  NY  14221 

1 Civil/Nuclear  Systems  Corp 
ATTN : Robert  Crawford 
1200  University  NE 
Albuquerque,  NM  87102 

1 EG§G,  Incorporated 
Albuquerque  Division 
ATTN:  Technical  Library 
P.  0.  Box  10218 
Albuquerque,  NM  87114 

1 General  American  Trans  Corp. 
General  American  Res  Div 
ATTN:  G.  L.  Neidhardt 
7449  N.  Natchez  Avenue 
Niles,  IL  60648 

1 General  Electric  Company-TEMPO 
ATTN:  DAS I AC 
P.  0.  Drawer  QQ 
Santa  Barbara,  CA  93102 


1 President 

General  Research  Corporation 
ATTN:  Library 

McLean,  VA  22101 

1 Honeywell,  Inc. 

Govt.  $ Aero  Products  Div 
600  Second  Street,  NE 
Hopkins,  MN  55343 

1 Hughes  Helicopters 
ATTN:  A.  C.  Edwards 
Dept.  15-43 

Centinela  5 Teale  Streets 
Culver  City,  CA  90230 

1 J.G.  Engineering  Research 
Associates 
3831  Menlo  Drive 
Baltimore,  MD  21215 

1 J.L.  Merritt  Consulting  § 

Special  Eng  Services,  Inc. 
ATTN:  Technical  Library 
P.  0.  Box  1206 
Redlands,  CA  92373 

2 Kaman  AviDyne 

ATTN:  Dr.  N.  P.  Hobbs 
Mr.  S.  Cricione 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01803 

2 Lockheed  Missiles  $ Space  Co. 
ATTN:  T.  Gerrs  D/52-53 
Bldg.  205 

Tech  Info  Ctr,  Doc/Coll. 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

1 McDonnell  Douglas  Corporation 
ATTN:  Robert  W.  Halprin 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
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1 The  Mitre  Corporation 
ATTN;  Library 

Rt  62  $ Middlesex  Turnpike 
P.  0.  Box  208 
Bedford,  MA  01730 

2 Physics  International  Corp. 
ATTN:  E.  T.  Moore 

Dennis  Orphas 
2700  Merced  Street 
San  Leandro,  CA  94577 

4 Physics  International  Corp. 
ATTN:  Robert  Swift 

Charles  Godfrey 
Larry  A.  Behrmann 
Technical  Library 
2700  Merced  Street 
San  Leandro,  CA  94577 

5 R§D  Associates 

ATTN;  Dr.  H.  L.  Brode 

Dr.  Albert  L.  Latter 
Bruce  Hartenbaum 
William  B.  Wright 
Henry  Cooper 
P.  0.  Box  9695 
Marina  del  Rey,  CA  90291 

4 R£jD  Associates 

ATTN:  Jerry  Carpenter 
Sheldon  Schuster 
J.  G.  Lewis 
Technical  Library 
P.  0.  Box  9695 
Marina  del  Rey,  CA  90291 

1 Sandia  Laboratories 

ATTN:  Doc  Control  for  3141 
Sandia  Report  Coll. 
Albuquerque,  NM  87115 


2 Sandia  Laboratories 
Livermore  Laboratory 

ATTN : Doc  Control  for  Tech  Lib 
Doc  Control  for  L.  Hill 
P.  0.  Box  969 
Livermore,  CA  94550 

1 Science  Applications,  Inc. 

ATTN:  R.  A.  Shunk 

P.  0.  Box  3507 
Albuquerque,  NM  87110 

3 Science  Applications,  Inc. 

ATTN:  R.  Seebaugh 

William  M.  Layson 
John  Mansfield 
1651  Old  Meadow  Road 
McLean,  VA  22101 

2 Science  Applications,  Inc. 

ATTN : David  Bernstein 

D.  E.  Maxwell 

2450  Washington  Ave,  Suite  120 
San  Leandro,  CA  94577 

2 Science  Applications,  Inc. 

ATTN:  Technical  Library 
Michael  McKay 
P.  0.  Box  2351 
La  Jolla,  CA  92038 

1 Shock  Hydrodynamics,  Inc. 

A Division  of  Whittaker  Corp. 
ATTN:  L.  Zernow 

4710-16  Vineland  Avenue 
North  Hollywood,  CA  91602 

5 Systems,  Science  5 Software 
ATTN:  Robert  T.  Allen 

Donald  R.  Grine 
Ted  Cherry 
Thomas  D.  Riney 
Technical  Library 
P.  0.  Box  1620 
La  Jolla,  CA  92037 
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2 Terra  Tek,  Inc. 

ATTN:  Sidney  Green 

Technical  Library 
420  Wakara  Way 
Salt  Lake  City,  UT  84108 

2 Tetra  Tech,  Inc. 

ATTN:  Li-San  Hwang 

Technical  Library 
630  North  Rosemead  Blvd. 
Pasadena,  CA  91107 

7 TRW  Systems  Group 

ATTN:  Paul  Lieberman 

Benjamin  Sussholtz 
Norm  Lipner 
William  Rowan 
Jack  Farrell 
Pravin  Bhutta 
Tech  Info  Ctr/S-1930 
One  Space  Park 
Redondo  Beach,  CA  92078 

2 TRW  Systems  Group 

ATTN:  Mr.  F.  A.  Pieper 
Greg  Hulcher 

San  Bernardino  Operations 
P.  0.  Box  1310 
San  Bernardino,  CA  92402 

2 Union  Carbide  Corporation 

Holifield  National  Laboratory 
ATTN:  Doc  Control  for  Tech  Lib 

Civil  Defense  Res  Proj 
P.  0.  Box  X 
Oak  Ridge,  TN  37830 

1 Universal  Analytics,  Inc. 

ATTN:  E.  I.  Field 
7740  W.  Manchester  Blvd. 

Playa  del  Rey,  CA  90291 

1 URS  Research  Company 

ATTN:  Technical  Library 

155  Bovet  Road 

San  Mateo,  CA  94402 
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1 Weidlinger  Assoc.  Consulting 
Engineers 

ATTN:  J.  Isenbert 

2710  San  Hill  Road 
Suite  230 

Menlo  Park,  CA  99025 

1 Westinghouse  Electric  Co. 
Marine  Division 

ATTN:  W.  A.  Volz 
Hendy  Avenue 
Sunnyvale,  CA  94008 

4 Denver  Research  Institute 
University  of  Denver 
ATTN:  Mr.  J.  Wisotski 
Fred  P.  Vanditti 
Ron  W . Buchanon 
Technical  Library 
P.  0.  Box  10127 
Denver,  CO  80210 

3 IIT  Research  Institute 
ATTN:  Milton  R.  Johnson 
R.  E.  Welch 
Technical  Library 
10  West  35th  Street 
Chicago,  IL  60616 

2 Lovelace  Foundation  for 

Medical  Education 
ATTN:  Asst  Dir  of  Research 
Robert  K.  Jones 
Technical  Library 
5200  Gibson  Blvd,  SE 
Albuquerque,  NM  87108 

1 Massachusetts  Institute  of 
Technology 

Aeroelastic  and  Structures 
Research  Laboratory 
ATTN:  Dr.  E.  A.  Witmer 
Cambridge,  MA  02139 
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2 Southwest  Research  Institute 
ATTN:  Dr.  W.  E.  Baker 
A.  B.  Wenzel 
8500  Culebra  Road 
San  Antonio,  TX  78206 

4 Stanford  Research  Institute 
ATTN:  Dr.  G.  R.  Abrahamson 

Carl  Peterson 
Burt  R.  Gasten 
SRI  Library,  Rm  G021 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


2 Washington  State  University 
Administration  Office 
ATTN:  Arthur  Miles  Hohorf 
George  Duval 
Pullman,  WA  99163 

1 University  of  Tennessee 
Space  Insitute 
ATTN:  Dr.  D.  R.  Keefer 
Tullahoma,  TN  37388 


Aberdeen  Proving  Ground 


1 University  of  Illinois 
Consulting  Eng.  Services 
ATTN:  Nathan  M.  Newmark 
1114  Civil  Engineering  Bldg. 
Urbana,  IL  61801 

2 The  University  of  New  Mexico 
The  Eric  H.  Wang  Civil 

Engineering  Rsch  Facility 
ATTN:  Larry  Bickle 
Neal  Baum 

University  Station,  Box  188 
Albuquerque,  NM  87131 


Dir,  USAMSAA 

ATTN:  Dr.  J.  Sperrazza 

Mr.  R.  Norman,  GWD 
Dr.  Rivello 
R.  Bailey 

DRXSY-MP , H.  Cohen 
Cdr , USATECOM 

ATTN:  DRSTE-SG-H 
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